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^7^; Abstract 

^ • We investigate a steady flow of a viscous compressible fluid with inflow boundary condi- 

■ tion on the density and inhomogeneous slip boundary conditions on the velocity in a cylindrical 
domain Q = Qq x (0, L) G M^. We show existence of a solution (v.p) G Wp{Q) x Wp{Q), 
where v is the velocity of the fluid and p is the density, that is a small perturbation of a constant 
^ I flow (v = [1, 0, 0], p = 1). We also show that this solution is unique in a class of small pertur- 
^ • bations of (v, p). The term u ■ Vw in the continuity equation makes it impossible to show the 
existence applying directly a fixed point method. Thus in order to show existence of the solution 
^ . we construct a sequence (i;", p") that is bounded in Wp{Q) x Wp{Q) and satisfies the Cauchy 



condition in a larger space Loo(0, L; L2(fio)) what enables us to deduce that the weak limit of a 

O ' subsequence of (v", p") is in fact a strong solution to our problem. 
ON ' 



1 Introduction 

The mathematical description of a flow of a viscous, compressible fluid usually lead to problems 
of mixed character as the momentum equation is elliptic (in stationary case) or parabolic (in 
case of time-dependent flow) in the velocity, while the continuity equation is hyperbolic in the 
density. Therefore, the application of standard methods usually applied to elliptic or hyperbolic 
problems fails in the mathematical analysis of the compressible flows and a combination of 
such techniques, as well as development of new mathematical tools is required. As a result a 
consistent theory of weak solutions to the Navier - Stokes equations for compressible fluids has 
been developed quite recently in the 90' s, mainly due to the work of Lions [QTI and Feireisl O. 
An overview of these results is given in the monograph lfT6l . A modification of this approach in 
case of steady flows with slip boundary conditions has been developed by Mucha and Pokomy 
in a dwo dimensional case in [14] and in 3D in f23\. 

The issue of regular solutions is less investigated and the problems are considered mainly with 
Dirichlet boundary conditions. If we assume that the velocity does not vanish on the boundary, 
the hyperbolicity of the continuity equation makes it necessary to prescribe the density on the part 
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of the boundary where the flow enters the domain. In [|25B Valli and Zajaczkowski investigate a 
time-dependent system with inflow boundary condition, obtaining also a result on existence of a 
solution to stationary problem. The existence of regular solutions to stationary problems with an 
inflow conditon on the density has been investigated by Kellogg and Kweon [8] and Kweon and 
Song [fTOl. Their results require some smallness assumptions on the data, and the regularity of 
solutions is a subject to some constraints on the geometry of the boundary near the points where 
the inflow and outlow parts of the boundary meet. In [9] Kellogg and Kweon consider a domain 
where the inflow and outflow parts of the boundary are separated, obtaining regular solutions. 

The lack of general existence results inhibits the development of qualitative analysis of com- 
pressible flows. Therefore it is worth to mention here the papers by Plotnikov and Sokolowski 
who has investigated shape optimization problems with inflow boundary condition in 2D [|2T1l 
and 3D ll22l dealing with weak solutions. More recently Plotnikov, Ruban and Sokolowski have 
investigated shape optimization problems working with strong solutions in [19] and [2^. 

It seems interesting both from the mathematical point of view and in the eye of applications to 
investigate problems with inflow boundary condition on the density combined with slip boundary 
conditions on the velocity, that enables to describe precisely the action between the fluid and the 
boundary. Such problem is investigated in this paper. The domain is a three dimensional cylinder 
and we assume that the fluid slips along the boundary with a given friction coefficient and there 
is no flow across the wall of the cylinder. We show existence of a regular solution that can be 
considered a small perturbation of a constant solution. The method of the proof is outlined in 
the next part of the introduction and now we are in a position to formulate our problem more 
precisely. 

The flow is described by the Navier-Stokes system supplied with the slip boundary conditions 
on the velocity. The complete system reads 

pv ■ Vf — pAv — (/i + u)Vdw V + Vn{p) = in fi, 

div (pv) = in fi, 

n ■T{v,tt{p)) ■ Tk + fv ■ Tk = bk, A; = 1,2 on T, (1.1) 

n ■ V = d on r, 

P = Pin on Tin, 

where w : — > is the unknown velocity field of the fluid and p : ^ M is the unknown 
density. We assume that the pressure is a function of the density of a class C^. Further, p and u 
are viscosity coefficients satisfying (p + 2v) > and / > is a friction coefficient. The domain 
i7 is a cylinder in M.^ of a form = VIq x (0, L) where fio G is a set with a boundary regular 
enough and L is a positive constant (see fig. [U). We want to show existence of a solution that 
can be considered a small perturbation of a constant flow (v, p) = ([1, 0, 0], 1). Thus we denote 
the subsets of the boundary V = d^l as F = Fj„ U Tout U Fq, where Fj„ = {x E T : v ■ n < 0}, 
^out = {x eT : V ■ n > 0} andTo = {x eT : V ■ n = 0}. 

By n we denote the outward unit normal to F and ri , T2 are the unit tangent vectors to F. Since 
the boundary has singularities at the junctions of Fj„ and Tout with Fq, for the boundary traces 
we will consider functional spaces that are algebraic sums of spaces defined on the boundary. 
More precisely for s,q eRwe shaU denote W^{T) := W^{Tin) + W^{Tout) + ^^(ro)- We 
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assume that b e Wp ^^^{T), p,„ G W^{Ti^) and d G Wp ^^^(T) are given functions and d = 
on Tq what means that Fq is an impermeable wall. 

For simplicity we consider the momentum equation with zero r.h.s., but our proofs work 
without any modification for the r.h.s. p F where F is small enough in Lp. 

We shall make here some remarks concerning notation. Since we will usually use the spaces 
of functions defined on ^l, we will skip in notation of the spaces, for example we will write 
L2 instead of ^2(^2). For the density we will use estimates in the space Loo(0, L; L2(i7o))- For 
simplicity we will denote this space by L^(L2). A constant dependent on the data that can be 
controlled, but not necessarily small, will be denoted by C, and E shall denote a constant that 
can be arbitrarily small provided that the data is small enough. 

In order to formulate our main result let us define a quantity Dq that measures how the 
boundary data b,d and pin differ from the values of, respectively, fv-Ti,n-v and p in appropriate 
norms. We have v ■ Ti = r^^^ and v ■ n = n^^^ thus we define 

^0 = \ h- fr-^'^Ww^-i/p^^r) + W^-^^^'^Ww^-'^^iF) + \ \Pin-M\w^(r,„)- (1-2) 
Our main result is 

Theorem 1. Assume that Dq defined in (1221) is small enough, f is large enough and p > 3. Then 
there exists a solution (f , p) G W^{Q) x W^iVL) to the system ( li.iP and 

+ \ \p-p\\w^ < E{Do), (1.3) 

where E{Dq) can be arbitrarily small provided that Dq is small enough. This solution in unique 
in the class of solutions satisfying the estimate U.3h 

The major difficulty in the proof of Theorem H] is in the term u ■ Vw in the continuity equa- 
tion, that yields impossible a direct application of a fixed point argument. To overcome this 
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problem one can apply the method of elliptic regularization, known rather from the theory of 
weak solutions (see [T6]). This method has been applied to a similar problem in a two dimen- 
sional case in [18]. However, it complicates considerably the computations since we have to find 
the bound on the artificial diffusive term. Here we apply a method of successive approximations, 
that leads to a more direct proof. In order to prove Theorem \T\ we will construct a sequence 
(m", w") G Wp X Wp that converges to the solution of (|l.ll) . Due to the presence of the term 
u ■ Vw we can not show directly the convergence in Wp x Wp, but we can show that (m", w"^) 
is a Cauchy sequence in a larger space x L^{L2) and thus converges in this space to the 
weak solution of (11.11) . On the other hand, the sequence will converge on a subsequence weakly 
in Wp X Wp, what will enable us to show that the weak solution is in fact strong. A similar ap- 
proach has been applied in [Q to an evolutionary Navier-Stokes system in a framework of Besov 
spaces. 

We start with removing the inhomogeneity from the boundary condition (ll.ll) /i. To this end 
let us construct mq G Wp{VL) such that 



n ■ uqW 



d-rS^\ (1.4) 
Due to the assumption of smallness of d — n^^^ in Wp~^^^(T) we can assume that 

\\uq\\w^ « 1. (1.5) 
From now on we assume (11.51) in all our results. Now we consider 

u = V — V — Uq and w = p — p. 
One can easily verify that (u, w) satisfies the following system: 

dx-^u — piAu — {u + p)'Vdivu + n'{l)\/w = F{u,w) in 

div u + dxj^w + {u + Uq) ■ Vw = G{u,w) in Q, 

n-2fiB{u) ■ Ti + f u- Ti = Bi, i = l,2 on T, (1.6) 

n - u = on r, 

w = Wi„, on Vi 



mi 



where 



F(m, w) = —W {u + V + Uq) ■ V(tt + Uq) — {Uq ■ Vu) — U ■ Vtto 

+pAuQ + (z/ + ^)Vdivwo - uq ■ Vuq - [tx'{w + 1) - 'k'{1)]Vw, (1.7) 
G{u, w) = —{w + 1) divuo — wdivu 



and 



Bi = bi-2fin- B{uo) ■ n - /r, 



(1) 



From now on we will denote 7r'(l) =: 7. We see that F and G also depend on Vm, uq, Vuq, 
but for simplicity we will write F{u, w) and G{u, w). In order to prove Theorem[T]it is enough 
to show the existence of a solution {u,w) to the system (11.61) provided that 1-1/p, , and 
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I l^ol |vK2(r2) are small enough. As we already mentioned, we will construct a sequence that con- 
verges to the solution. The sequence will be defined as 

- /iAu"+i - (z/ + /i) Vdiv + 7Vu;"+^ = F(u", w"") in Q, 

divM"+i + d^,w''+^ + (u" + uo) ■ Vw''+^ = G{u'\ w'') in fi, 

n ■ 2/iD(M"+i) ■Ti + f ■ Ti = i = 1, 2 on T, (1.8) 

n • u'^+i = on r, 

= y;.^ On Tin- 

As we will see in the sequel, our method does not require any particular starting point for the 
sequence but only some smallness assumptions on the starting point hence 

without loss of generality we can set w^) = (0, 0). In order to show the existence of the 
sequence defined in (11.81) we have to solve a linear system: 

dx-t^u — fiAu — (z/ + /i) Vdiv u + 'yVw = F in 

div u + dx^w + {u + uq) ■ Vw = G in Q, 

n-2nT){u)-Ti + fu-Ti = Bi, i = l,2 on T, (1.9) 

n ■ u = on r, 

W = Win on Tin, 

where (F, G, u, uq) E Lp x x Wp x Wp are given functions and ■ n = on F. 

Let us now outline the strategy of the proof, and thus the structure of the paper. In section [2] 
we show the a priori estimate (12.351) on a solution to the linear system (11.91) . We start with an 
energy estimate in x Loo(i^2)- Next the properties of the slip boundary conditions enables us 
to show that the vorticity of the velocity on the boundary has the same regularity as the velocity, 
and this fact makes it possible to find a bound on | Then the estimate (12.351) results 
directly from the elliptic regularity of the Lame system. 

The linear system (11.91) is solved in section [3l First we show the existence of a weak solution 
using the Galerkin method modified to deal with the continuity equation. Next we can show 
that this solution is in fact strong using a priori estimate and symmetry of the slip boundary 
conditions. 

In section|4]we show the estimate in x Wp on the sequence (w*^, ly") and, as a result, the 
Cauchy condition satisfied by this sequence in the space x Loo(i^2)- These results are derived 
by application of the estimates for the linear system. 

In section [5] we apply the results of section |4] passing to the limit with {u"',w''^) and then 
showing that the limit is a solution to (11.61) . Finally we show that this solution is unique in a class 
of solutions satisfying the estimate (11.31) . 



2 A priori bounds 

The main result of this section is the estimate (12.351) in Wp x Wp. In order to show it we start 
with an energy estimate in x Lca{L2). Next we consider the equation on the vorticity of the 
velocity and apply the Helmholtz decomposition to derive the bound on | |iy | and finally using 
the classical elliptic theory we conclude (12.351) . 
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In our proofs we shall not need explicit formulas on the functions F{u, w) and G{u, w), what 
will be important is that they depend quadratically on u and w. More precisely, we will show a 
following estimate 

Lemma 1. Let {u, w) G Wp x and let F{u, w) and G{u, w) be defined in f li.TD . Then 

\\F{uM\W + \\G{u,w)\\wl< .2 1) 

Proof Since by the imbedding theorem (^^ ) C Loo (^) , the estimate on 1 1 G 1 1 14/^1 is straight- 
forward, and the only part of F that deserves attention is 5tx\w)Ww, where 

Jtt'H :=7r'(^i; + l)-7r'(l). (2.2) 

We will apply a fact that for a - function / we have 

/(x) - f{y) = {x-y) C f[tx + (1 - t)y] dt, (2.3) 

Jo 

Thus we have 

»i 



6n'{w) = w n"{tw + 1) dt. 
Jo 

Since tt is a - function, the above implies 

\\6n'{w)Vw\\L, < C{7r)\\w\U\Vw\\L, <C\\w\\l.r. 

The other parts of F can be estimated direcly giving (|2.1I) . □ 

Next, we derive the 'energy' estimate in x ^00(^2). It is stated in the following lemma 

Lemma 2. Let {u,w) be a solution to the system di.9D with {F,G, B,Win,u) & V* x L2 x 
L2(T) X L2{Tin) X Wp, with \ \u\\w^ small enough and f large enough. Then 

M\m + \Ml^{l,) <C[\\F\\v* + \\G\\l, + \\B\\L,ir) + IkmlU^cr^J, (2.4) 

where 

V = {v e H\n) -.vnlr = 0} (2.5) 

and V* is the dual space ofV. 

Proof. We apply a general identity 



J^{—fiAu — {1/ + /i)Vdivu) vdx = 
= J^{2yuD(w) : V f + z/div u divf } dx — J-^n ■ [2fi'D{u) + u div -uld] ■ v da. 

For u, V satisfying the boundary conditions (11.91 ) c; 4 the boundary term in (12.61) equals 

[{f2[B^-fiu■n)]iv■n)}da. 



(2.6) 



Thus multiplying (|1.9l) i by u and integrating over we get 



J^{2 fiD"^ (u) + udiv'^ u} dx + Jj.{f + ^)\u\'^ da — j J^w div udx = ^^q-^ 
= j^F ■ udx + J^{Bi{u ■ Ti) + B2{u ■ T2)} da. 

From now on (not only in this proof but also later) we will use the summation convention when 
taking the sum over the tangential components. Applying (11.91 )9 and the boundary conditions we 
get 

/ wdivudx= / Gw dx + - / w"^ div {u + uq) dx 
Jn Jn 2 Jq 

wHl + u'i^)da+l- [ wl{l + u'i^)da. 



2 Jr . 2 

out ^ ■ 

For ||uo||vy2 small enough we have by the imbedding theorem 1 + Ug^-* > a.e. on V out what 

yields Jp 10^(1 + da > 0. Moreover, for the friction / large enough on Tin the boundary 
term in (12.71) will be positive. Combining these facts with the Kom inequality (that can be proved 
in a simple way with the friction large enough - see Lemma 2.4 in IfTTIl ): 

[ 2fi-D\u)+ I f{u-rfda>C\\u\\]j, (2.8) 

we derive from (12.71) the following inequality 

C||it||/i-i < / F-udx+ I Bi{u-Ti) da+- f w"^ div (u+uq) dx—- f wfni^+u'i^) da. (2.9) 
Jn Jr ^ Jn 2 Jp.^ 

In order to derive (12.41) from (12.91 ) we have to estimate | |loo(L2) i^i terms of and the 

data. To show this estimate we refer to section [3] where the linear system (11.91) is solved. Namely, 
we have w = S{G — div u) where the operator S is defined in (|3.7I) and thus the estimate (13.81) 
implies 

\\w\\l^{L2) <C{\\G\\l2 + \\u\\h^ + lkm||L2(r.„))- (2.10) 

The above inequality combined with (12.91 ) yields (12.41 ). □ 

Now we consider the vorticity of the velocity a = rot u. The properties of the slip boundary 
conditions enables us to express the tangential components of a on the boundary in terms of the 
velocity. We arrive at the following system 



dx^a — fiAa = rot F in fi, 

a-T2 = {2xi~i)u-n + f on T, 

a.n = {l-2x2)u-T2-f on T, ^^'''^ 

div a = on r, 



where Xi denote the curvatures of the curves generated by tangent vectors r^. In order to show 
the boundary relations (12.1 11 )9 it is enough to differentiate (11.91 )4 with respect to the tangential 
directions and apply (II. 9h . A rigorous proof, modifying the proof in the two-dimentional case 
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from [[TSl . is given in the Appendix. The condition div a = in i7 results simply from the fact 
that a = rot u. We introduce this relation as a boundary condition (|2.11I) a. that completes the 
conditions on the tangential parts of the vorticity. What is remarkable in the boundary conditions 
(12.1 11) 9;:! is that the tangential parts of the vorticity on the boundary has the same regularity as the 
velocity itself and the data. This feature of slip boundary conditions makes it possible to show 
the higher estimate on the vorticity (see [[l2l,[[n]|, [l23ll '). 

In order to derive the bound on the vorticity we can follow (23), Lemma 4, and construct 
ao, a divergence-free extension of the boundary data fl2.1 ip o for example as a solution to the 
Stokes problem with zero r.h.s and the boundary conditions (12.1 1|) 9 -:! supplied with ao ■ n = 0. 
The theory of the Stokes system then yields 

||«o||H/pi ^ C [\\u\\^i-i/p^^^ + \ \B\\^i-i/p^^^]. (2.12) 

Then the function a — ao satisfies the system 

—fiA{a — ao) = rot [F — dx^u] + fiAao in Q, 

{a - ao) ■Ti = on T, 

{a — ao) ■ r2 = on F, 

div {a — ao) = on F. 

Here we have used the fact that dx^a = lotdx^u to preserve the rotational structure of the r.h.s. 
For the above system we have the following estimate (see [|26l ') 

< C [\\F\\L, + \\dxM\L, + \\ao\\w^]- (2.14) 

The term with ao can be bounded by (|2.12l) and to deal with dxj^u we apply the interpolation 
inequality (|6.3I) . We obtain the term | l^i that we bound using (|2.4I) and finally arrive at 

\\a\\w^ < C(e) [\\F\\lp + \ \G\\wi + lh»n||L2{r,„) + 1 1^1 + 1 15| |^i-i/p^j.^] + e||u||i^^2. 

(2.15) 

With the bound on the vorticity at hand the next step is to consider the Helmholtz decomposition 
of the velocity (the proof can be found in [7 |): 

u = V(j) + A, (2.16) 

where (f)\r = and div A = 0. We see that the field A satisfies the following system 

rot A = a in Vl, 

div^ = in fi, (2.17) 
A-n = on F. 

This is the standard rot-div system and we have ||74||iy2 < C | |q!| what by (|2.15l) can be 
rewritten as 

\\A\\w^ < C(e) [\\F\\l^ + \ \G\\wi + \ \u\\^i-i/p^^^ + | |S| + \ \win\\w}(r,J + e\\u\\w^ 

(2.18) 
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for any e > 0. Now we substitute the Helmholtz decomposition to (|1.9l )i . We get 

V[-{u + 2^)A(/) + -fw] = F ~ d^,A + nAA + {u + fi)V div A - d^,(f}, (2.19) 
but = div u and denoting the l.h.s. of the above equation by F we obtain 

+ 2^)div u + -fw = H, (2.20) 
where VH = F. Combining the last equation with (11.91) 9 we arrive at 

7W + + {u + uq)Vw = H, (2.21) 

where 7 = —V and 

H = — ^ + G. (2.22) 

The equation (12.211) makes it possible to estimate the VF^ -norm of the density in terms of Wp - 
norm of H. The latter will be controlled since (12.191) enables us to bound 1 1 Vi/| II^ and | | |lp 
using interpolation and the energy estimate (12.41) . The details are presented in the proof of lemma 
HI but first we estimate | | ^4/^1 in terms of H. The result is stated in the following lemma 

Lemma 3. Assume that w satisfies the equation ([2.21\) with H G Wp. Then 

\\w\\w:^ < C [\\H\\w:^ + |km|k^i(r,„)] • (2.23) 

Proof. In order to find a bound on | |if | II^ we multiply (12.211 ) by \w\^~'^w and integrate over 
VL. Integrating by parts and next using the boundary conditions we get 

I \w\^~'^w Wx^dx = - I dx-^\w\^ dx = - I \wf da f \w\^ da, 

Jn P Jn P Jvout P JTin 

since n*^^) = on Fq, n^^^ = —1 on Tin and n^^^ = 1 on Tout- Similarily, applying the boundary 
conditions we get 

I {u + Uq) ■ {\w\^~'^w'Vw) dx = - I {u + uo) -Vlwl^ dx = 
Jn P Jn 

1 ^ - ■ 1 ^ ,« L.,iP^^ 1 f „,a) 



(2.24) 



div {u + uq) \w\^ dx -\ — / Uq \w\^ da / Uq \w\^ da. 

P Jn P Jvout P Jvi^ 

Thus multiplying (12.211) by \ w\^~'^'w we get 

7l hi II, - I In (u + uo) \w\p dx + i /^^^^ \w\p (1 + 4^^) da < 
<\\H\\LjHrC + lJ,Jw,n\Hl + nl^^)da. 

By the imbedding theorem the smallness of | |m + mo| Ivk^ implies 1 + Mq^^ > a.e. in Q and 
7 — ||div (u + Mo) I loo > 0. Thus the boundary term on the l.h.s. is positive and the term with 
div {u + uq) can be combined with the first term of the l.h.s, what yields 

Clltyllr < Ili^fllL I |ty| I r""^ + C I IWjnl 1^ - 



and so 

lkllL,<C[||//|U, + IK.|lMr.„)]- (2.25) 

The derivatives of the density are estimated in a similar way. In order to find a bound on Wx^ we 
differentiate (12.211) with respect to Xi. If we assume that w G then (12.211) implies u ■ Vw G 
Wp, where 

u:=[l + {u + UoY^\ {u + UoY^\ {u + uoY% (2.26) 

Thus u ■ Vwxi ■= (u ■ Vw)x^ — ■ Vw G Lp. Hence we can differentiate (12.211) with respect to 
Xi, multiply by \wxi\^~'^Wxi and integrate. Since Ux^ = {u + uo)xi, we have 

/ Uxr {\wxf~^Wx,Vw) dx < \\V{u + uo)\\l^ ||Vt/;||^ < C \\u + uo\\w^\\Vw\\l,. 

Next, since w ■ Vwxi G Lp, we can write 

/ u-\wx\^~'^WxVwxdx=- [ u-VliVxl^ dx = — / \wx\^ diYudx+- [ \wx\^u-nda = 
Jn pJn ' pJn ' P Jr ' 

= I \wx^\^ divudx f \win,x,\^{l + u'(^^)da+- f |w^J^ (1 + 4^^) rfa. 
P Jn P JTi^ P JTout 

For z = 2, 3 we have Wj„,x, ^ LpiVin) and hence the above defines the trace of \wxS^ on ^out- 
We arrive at 



\w, 



(2 27) 



For i = 2, 3 it gives directly the bound on \ \wx^ \\l ■ In order to estimate Wx^ we start the same 
way differentiating (12.211) with respect to xi and multiplying by \wxi I'^^'^Wxj^. The difference in 
comparison to Wx^ and Wx^ is that Wx^ is not given on Tin. In order to overcome this difficulty 
we can observe that on Tin the equation (12.211) reduces to 

JWin + {U + UoY^^ Win,x2 + {u + UqY^^ Win,x3 + [! + (« + UqY^^] Wx^ = H, 

what can be rewritten as 

1 



1 + (m + mo)« 

Thus we have 



[H - -fWin - {U + Uo)r " Vr^m] • 



IkxilUpCr^) < C [\\H\r^J\Lp(r,„) + \ \win\\w^{r,„)]- 
Using this bound in (12.271) . i = 1, we arrive at the estimate 

< C [\\Hx,\\l, \M\l~^+\\u+uo\\w2\\Vw\\l^+\\H\\l^^^^^^^ (2.28) 
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The boundary term | |Lp(ri„) can by replaced by due to the trace theorem. Thus 

combining (12.271) (for x-z and X3) with (|2.28l) we get 

\\Ww\\1^<C[\\VH\\l,\\Ww\\1-' + \\u + uo\\wi\\Vw (2.29) 



The term | |ii + iio| liya | \Vw\ can be put on the l.h.s. due to the smallness assumption and thus 
we get 

\\^w\\l, <C[\\H\\w^^ + |h.„||H/,i(r,„)], (2.30) 
what combined with (|2.25l) yields 

\\w\\w^ < C [\\H\\wi + \\H\\L^^r^„) + \\win\\w^ir,„)]- (2.31) 

Applying again the trace theorem to the term \ \H\ |Lp(r,„) we arrive at (12.231) . □ 
The next step is to estimate H in terms of the data. The result is in the following 

Lemma 4. Let H be defined in i[2.22\l . Then ^6 > we have 

\\H\\w^i < 5\\u\\w^ + C{5)[\\F\\l^ + \\G\\w} + \\B\\^^^^/,^^^ + \\wJ\wi(r,J. (2.32) 
Proof. Applying first the interpolation inequality (16.31) and then the estimate (12.41) we get 

Ili^lU, < 5i||V//|U, + Ci6^) [\\F\\l, + \\G\\l, + \\B\\L,ir)]. (2.33) 
Next, by (12.191) we have 

llVi^llr < C [\\F\\r + ||G||wi + ||A||w. + WO^MlI 



where u = Vcf) + A is the Helmholtz decomposition. Now we use the bound (12.181) on 1 1 |vy2. 
We obtain a term 1-1/p, that we estimate using the trace theorem and the interpolation 
inequality (16.31 ). The same inequality is applied to estimate | l^^^j^l l^.^. We arrive at 

\\VH\\l^ < C [\\F\\l^ + \ \G\\wi + , + ||wm||vyi(r>n)] 



+6,\\u\\wi + C{5,)[\\F\\l, + \\G\\l, + II^IImd]. 



(2.34) 



Combining (12331) and (12341) we get (12321) □. 

Now we are ready to show the a priori estimate in x on the solution of the linear 
problem. 

Lemma 5. Let {u, w) be a solution to ( li.9D with {F, G, B, Win, u) e Lp x x wl'^^^iV) x 
Wp{Tin) X Wp, with I |m| I small enough and f large enough. Then 

M\w^ + \H\w^ <G[\\F\\l^ + \\G\\w^ + ||5||^^i-i/p(^^ + \\win\\w^(r,J- (235) 
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Proof. If (w, w) is a solution to (11.91 ). then in particular the velocity satisfies the Lame system 

dx^u — fiAu ~ {u + fi)'Vdiv u = F — '-f'Vw in Q, 

n ■ 2^D(u) ■Ti + f u-Ti = Bi, 1 = 1,2 on T, (2.36) 
n ■ u = on r. 

The classical theory of elliptic equations (Agmon,Douglis,Nirenberg ll2l,[[3]|) yields 

\\u\\w^ < C[\\F\\l^ + \ \w\\wi + + \ \u\\wi]- 

Applying the interpolation inequality (16.31) to the term | |m| li^/^; and then the energy estimate (12.41) 
we get 

\\u\\w^ < C[||F||l^ + \ \G\\w} + \ \w\\wi + + ||wi„||L2(r,„)]. (2.37) 

In order to complete the proof we combine (12.231) and (12.321) obtaining 

\\w\\w^ < 6\\u\\w^ + C{6)[\\F\\l^ + \\G\\w^ + + \\win\\w}(r,J, (2.38) 

and choosing for example ^ = ^ where C is the constant from (12.371) we arrive at (|2.35l) . □ 

3 Solution of the linear system 

In this section we show the existence of the sequence (u", w"') defined in (11.81) . To this end we 
have to solve the linear system (11.91) where {F, G, u, uq) & Lp x x x are given 
functions such that w ■ n = on F. First we apply the Galerkin method to prove the existence 
of a weak solution and next we show that this solution is strong. For simplicity we will denote 

M + Wo by u. 

3.1 Weak solution 

Let us recall the definition of the space V (12.51 ). A natural definition of a weak solution to the 
system (11.91 ) is a couple [u, w) eV x Loo(i^2) such that 



{v ■ dxiU + 2^T){u) \W V + udwudwv — ^wdiYv} dx + I f{u-Ti){v-Ti)da = 
n Jr 

= / F-vdx+ / Bi{vr,)da (3.\) 
Jq Jr 

is satisfied W v E V and ([111)2 is satisfied in V'{n), i.e. W (p E C°°{ny. 

— / wu-'Vcpdx— / W(pdivudx+ / w(f)da= (l){G — div u) dx + / Win4>da, (3.2) 
Jq. Jn Jvout Jo Jrir, 
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where u is defined in (12.261) . Let us introduce an orthonormal basis of V: {uji]°^^. We consider 
finite dimensional spaces: = {^^^ ctiUJi : G R} C V. The sequence of approximations 
to the velocity will be searched for in a standard way as = X^ili ^i- Due to the equation 
(11.91) 9 we have to define the approximations to the density in an appropriate way. Namely, we set 
^ g(^QN _ divu^), where S : L2(f^) ^ 1.00(^2) is defined as 

( dr,,w + u-Vw = v in V'{n), 
w = S{v) <^ ' ' ' (3.3) 

{ w = Win on Tj^. 

We want the image of S to be in the space Loo(-Z^2) so that we can apply the theory of transport 
equation treating xi as a 'time' variable to show that S is well defined. In order to solve the 
system on the r.h.s. of (13.31) we can search for a change of variables x = ■ip(z) satisfying the 
identity 

d,,=d^,+u-V^. (3.4) 
We construct the mapping ^ in the following 



Lemma 6. Let \ \u\\y^2 be small enough. Then there exists a set U (ZM? and a dijfeomorphism 
X = il){z) defined on U such that Vt = '0(f/) and l\3.4\) holds. Moreover, if Zn ^ z and ilj{zn) — ^ 
To then n^{z) = 0, where n is the outward normal to U. 

Before we start with the proof we shall make one remark. The last condition states that the 
first component of the normal to ^^^(Fo) vanishes, but since ip is defined only on U we formulate 
this condition using the limits. It means simply that the image U = ip~^{il) is also a cylinder 
with a fiat wall. It will be important in the construction of the operator S. 

Proof of lemma m The identity (13.41) means that must satisfy 

^ = ^ = f.V). ^ = (3.5) 

OZi OZi OZi 

A natural condition is that V^(rj„) = rj„. Thus we can search for V'l-Zii -2^2, 2:3) = iIj^^^zA^i), 
where for all (2:2, 2:3) such that (22, -23, 0) G Tin the function ■022,^3 (') is a solution to a system of 
ODE: 



^'22,23(0) = (0,^2,^3)- 



(3.6) 



The r.h.s of the system (13.61) is a Lipschitz function with a constant K = ||Vm||oo and thus 
provided that ||m||h/2 is small enough the system (13.61) has a unique solution defined on some 
interval (0, 6^1,22): where b^^^^^ depends on Z2, z^ and ||Vm||oo- Provided that the latter is small 
enough the function ip(z) = ipz2,zz{zi) will be defined on U such that Vt = tp{U). 

Now we show that 'ip{z) = ^^22,23(21) ^ diffeomorphizm. The derivatives with respect to 
Zi are given by (13.51 ) and the remaining derivatives can be expressed in terms of u so we can see 
that J ijj = 1 + E{u), where E{u) is small (and thus > 0) provided that | |m| Ivv^ is small. 

To see that ^ is 1 — 1 we can write it in a form il){z) = z + e{z), where 1 1 Ve| is small. 
Assume that ip{z^) = ipiz"^) and z^ 7^ z"^. Then there exists i such that \zl — zf\ > \\z^ — z^] (the 



13 



lowercase denotes the coordinate). On the other hand, we have \zl — zf \ = \ei{z^) — £4(2^) | < 
1 1 1 1 Loo I ~ -2^^ I' what contradicts the smallness of 1 1 Ve| 

We have shown that the mapping ^jJ given by (|3.6I) is a diffeomorphizm defined on U such 
that ipiU) = ^l. Let us denote (p = tp^'^. Now it is natural to define the subsets of du as 
du = Uin U Uout U Uo where Uin = Tin, Uout = {z : z = lim(p{xn), Xn Tout} and 
Uq = {z : z = Um0(x„), x„ Tq}. 

In order to complete the proof we have to show that n^{z) = for 2; G Uq. But to this end it 
is enough to observe that 

Dij{z) {[1,0,0]) = [l+u\x),u\x),u'^{x)], 

where x = ipiz). But for x E Tq the vector on the r.h.s is tangent to Tq since u ■ n = 0. We can 
conclude that on Uq the image in of a straight line {{s, Z2, Z3) : s G (0, 6)} is a curve tangent 
to To, and thus Uq is a sum of such lines and so we have n^{z) = 0. The proof of lemma [6] is 
completed. □. 

Now we can define S{v) for a continuous function v as 

|-</>l(x) 

S{v){x) = Win{0,(f)2{x),(p3{x)) + v{i!{s , (p2{x) , 4>3{x))) ds . (3.7) 

Jo 

The condition = on 0(ro) guarantees that a straight line {s, zi, Z2) : s E (0, b) has a picture 
in and thus we integrate along a curve contained in ^l. It means that S is well defined for 
continuous functions defined on Q and the construction of ^ clearly ensures that S satisfies (|3.3I) . 
Next we have to extend 5" on L2(fi). To this end we need an estimate in Loo{L2)- It is given by 
the following 

Lemma 7. Let S be defined in ^3. 7\i . Then 

\\S{v)\\l^(L2) < C[\\Win\\L2iV,r.) + lkllL2(n)]. (3.8) 

Proof. Let Vl^^ denote an xi - cut of ^2 and let x := (x2, 2:3). Then by (|3.7I) we have 

\\S{v)\\\^^^^ ) = j [Win{0,<p2{x),<p3{x))+ j V{lp{s,(t)2{x),(t)3{x)))ds\ dx 

Jn^^ Jo 
< 2|km||i2(r,„) +C' / / v'^{Hs,Mx),Mx)))dsdx <C[\\win\\l^(^r,r^) + M\h{n)]^ 

J^x^ Jo 

The above holds for every xi E (0, L) what implies (|3.8I) . □ 

Now we can define S{v) for v E L2{VL) using a standard density argument. Let us take a 
sequence of smooth functions Vn ^ v'm L2{Vl). By (13.81) the sequence S{vn) satisfies 

\\S{vn)\\L^(L2) < C'[lkmlU2(r,„) + sup„||^;„||l2]. (3.9) 
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The bound on the r.h.s. is uniform in n and thus S{vn) ^* t] in Loo(i^2), and r] satisfies the 
estimate (|3.8I) . In particular for <j) E C°°{ri) we have 



S{vn)u ■ V(p dx ^ / rju-V(f)dx and / S {vn)(f) div udx ^ / ricjxiivudx. 
Jn Jq Jn 

In order to show that r] = S{v), i.e. r] solves the system on the r.h.s. of (13.31) we have to show 
that L S(vn) (bda ^ L neb da. To this end notice that the proof of lemma |7] implies in 

J- out V / ' J I 

particular that | | |L2(ro„t) satisfies the estimate (13.91) . Thus S{vn) ^ (in L2{Tout) for some 
C G L2{Tout), and in particular ^ S{vn)4> da ^ (cp da. We have to verify that r]\r^^t = (. 

This would not be obvious if we only had S{vn) E Loo{L2), but indeed the proof of lemma |7] 
implies a stronger condition that supremum (not only the essential supremum) of | \S{vn) \ \L2in^-^) 
is bounded, thus we must have ( = 77|ro„f We have shown that u ■ Vr] = v in V'{il), thus indeed 

T] = S{v). 

Having the operator S well defined we are ready to proceed with the Galerkin method. Taking 
F = F^,u = u^ = ^.cfui,v = ujk, k = l...Nandw = w^ = S{G^ - divu^) in (l3ll) . 
where and are orthogonal projections of F and G on V^, we arrive at a system of 
equations 

B^{u^,uJk) = 0, k = l...N, (3.10) 
where : V'^ V'^ is defined as 

B^{C^,v^) = J^{C^d,,v^ + 2fi-D{^^) : Vv^ + diY^^diYV^}dx 
-7 SiG"" - div e"^) div v"" dx + J^lf (e^^ -Tj)- Bi] {v^ ■ Tj) da - j^F^ ■ v"" dx. 

(3.11) 

Now, if satisfies (ITTOl) for A; = 1 . . . and ^x;^ = S{G^ - div u^), then a pair [u^ , w^) 
satisfies (13.11) - (13.21) for (f , 0) G (V^ x C'°°(fi)). We will call such a pair an approximate 
solution to (I37T]) - (IHl) . 

The following lemma gives existence of solution to the system (13.101) : 

Lemma 8. Let F,G E L'^{^), Win E L2(Tin), B E L2iT) and assume that f is large enough and 
\\u\\]y2 is small enough. Then there exists E satisfying ( li.iOP for k = 1 . . . N. Moreover, 

\\u^\\m <C{DATA). (3.12) 

Proof. In order to solve the system (|3.10l) we will apply a well-known result in finite- 
dimensional Hilbert spaces, lemma [14] in the Appendix. Thus we define the operator : 

_ as 

pN^^N^ = J2BN^^N^^^^^^ for e^GF^. (3.13) 

k 

In order to apply lemma [141 we have to show that {P{C^) , > on some sphere in . Since 
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B^(-, ■) is linear with respect to the second variable, we clearly have 

Jn Jq 



h 



Jn Jr Jn 



(3.14) 



h h 



j^F.edx-j^B,{e-n)da. 
Using the Korn inequality similarily as in the proof of the energy estimate (12.41 ) we get 

h + h>C\\e\\l, (3.15) 
for / large enough. We have to find a bound on /s. Denoting r]^ = S{G^ — div^^) we have 

- / rj^divC^dx= [ r]^{d,,r]^ + u-Vv^)dx- [ r]^ dx. (3.16) 
Jo Jn Jn 

Using (13.81) we get 

- fv''G'^dx> -llr/^IU. IIG^II^. > -CIIG^II^, + 11^^!^. + lk.n||L.(r.„)). 

Jn 

(3.17) 

With the first integral on the r.h.s of (|3.16l) we have 

In V^id^iV^ + u ■ Vt]^) dx = r]^{z)d,,ri^{z)Jij{z) dz = 

= J^l^^iz)^,,r^^iz) dz + J^^^iz)^,,r^^iz)[J^Piz) - 1] dz. ^^"'^^ 

The first integral can be rewritten as a boundary integral and since (z) = on (P{Tq), it reduces 
to 

If [r^'' {z)fn\z)da{z) = -\ j [r^'' {z)f da{z) + \ j [v"" (z)]' da{z) > - [ wlda{x). 

In the last passage we used the fact that (p\ri„ is the identity and that n^{z) > on Uouu what is 
true provided that cp does not differ too much from the identity on Touu what in turn holds under 
the smallness assumptions on u. 

With the second integral on the r.h.s. of (13.181 ) we have 



ri'\z)d,y {z)[J^{z) -l]dz> -sup^l - 1| / ii'\z){G'' - div^^ )(^) dz > 
u Ju 

> -i?lk''llL.(c/)[||G^||L.(f/)+||div,.e^|U,(^)] > -i?[||G'^||L(o)+ll?"^ll^i(o) + lh»|lL(r.„)]- 
Combining this estimate with (|3.15l) we get 

{P''{e),n > C m^'Wl.^n) - DWeWnnn) - D^], (3.19) 
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where D = \ \F\\L2^^) + \\G\\L2(n) + \\win\\L2{r„,) + \\B\\L2{r)- Thus there exists C = C{fi,Q,D) 
such that > for ||^|| = C, and applying lemma[l4]we conclude that 3^* : 

P^(^*) = and I iri I < C. Moreover, since {ujk}k=i is the basis of V^, we have P^(r) = 
(B^Ci^k) = 0, k = l...N. Thus C is a solution to (ITTOl) . □. 

Now showing the existence of the weak solution is straightforward. The result is in the 
following 

Lemma 9. Assume that F,G & L2{^1), Win ^ L2{Tin), B G -^2(r)- Assume further that f is 
large enough and \ [ yyg is small enough. Then there exists {u,w) G V x W that is a weak 
solution to the system 47. 91) . Moreover, the weak solution satisfies the estimate l \2.4\i . 

Proof. The estimates (D and (13121) imply that \ \u^\\m + \ \w^\\l^{L2) < C{DATA). 
Thus 

^ M in and ^* w in Loo{L2) 

for some (w, w) E x Loo{L2). It is very easy to verify that (tt, w) is a weak solution. First, 
passing to the limit in (13.11) for (m^, to^) we see that u satisfies (13.11) with w. On the other hand, 
taking the limit in (13.21) we verify that w = S{G — div u). We conclude that (m, w) satisfies (13.11) 
- (|3.2I) . thus we have the weak solution. To show the boundary condition on the density we can 
rewrite the r.h.s of (13.31) as 



+ tSiT^X2 + tS)^-3 = ttIw ^'(^)' (3.20) 
w = Win on r^, 

and, treating xi as a 'time' variable, adapt Di Pema - Lions theory of transport equation (f5]) that 
implies the uniqueness of solution to (13.201) in the class L^{L2). The proof is thus complete. □ 



3.2 Strong solution 

Having the weak solution of the linear system (11.91) we can show quite easily that this solution is 
strong if the data has the appropriate regularity. The following lemma gives existence of a strong 
solution to (fTgl) . 

Lemma 10. Let F e Lp, G e W^, Win G W^{Tin), B G Wp'^^^iT) and assume that f is large 
enough and \ \u\\w2 is small enough. Then there exist {u, w) G Wp x Wp that is a strong solution 
to ( |i.9l) and satisfies the estimate ([2.35\) . 

Proof. Since (11.91 ) is a linear system, the a priori estimate (12.351 ) will imply the regularity of 
the weak solution once we can deal with the singularity of the boundary at the juctions of Fq with 
Tin and Tout- This however can be done easily since Vl is symmetric w.r.t. the plane {xi = 0} and 
the slip boundary conditions preserve this symmetry. More precisely, for {x = (—xi,X2,Xs) : 
X = {xi,X2, Xs) G f2} we can can consider a vector field 

= [-u\x),u^{x),u^{x)]. (3.21) 

Then on Fj„ we have u ■ n = u ■ n and n ■ T){u) ■ ti + u ■ Ti = n ■ Ty{u) ■ ti + u ■ Ti. Hence we 
can extend the weak solution on the negative values of xi using (13.211) and, applying the estimate 
(12.351) , show that the extended solution is in Wp x Wp. An identical argument can be applied on 
Tout and we coclude that (w, w) is a strong solution to (11.91) . □ 
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4 Bounds on the approximating sequence 

In this section we will show the bounds on the sequence {(m", w")} of solutions to (11.81) . The 
term u ■ Vw in the continuity equation makes it impossible to show directly the convergence in 
Wp X Wp to the strong solution of (11.61) . We can show however that the sequence of iterated 
solutions is bounded in x W^, and using this bound we can conclude it is a Cauchy sequence 
in X Loo{L2), and thus converges in this space to some couple (u, w). On the other hand, 
the boundedness implies weak convergence in x W^, and the limit must be {u,w). The 
following lemma gives the boundedness of (u", to") in x W^. 

Lemma 11. Let {{u^, to")} be a sequence of solutions to UM starting from {vP , w^) = (0, 0). 
Then 

||«"lk,2 + lh1k,i<M, (4.1) 

where M can be arbitrarily small provided that ||uo||iy| (extension of the boundary data 
not to be confused with v!^ from [vP , w^), the starting point of the sequence w"^)), | |S| |^i 
1 1 '"'ml |wpi{r,„) <^nd I l^^l Iw"! <^re small enough and f is large enough. 

Proof. The estimate (|2.35l) for the iterated system reads 

|L,n+l|| „ I |L,,n+l|| ^ 

\\u I + 1 1^ I liyi S 
< C [\\Fiu-,w-)\\L, + \\Giu-,w-)\\w^. + ||S|Vi-i/.(r) + \M\w,Hr^„)]- ^"^'^^ 

Denoting A„ = | 1^/2 + | Iw-pi and 6 = \ \uo\\w2 + \\B\\^i-i/p^^^ + \\win\\w^(r,^),irom^J^ 
and (14.21) we get 

An+i <Al + b, (4.3) 

thus An is bounded by a constant that can be arbitrarily small provided that Aq and b are small 
enough. Indeed let us fix < 5 < | and assume that b < 6. Then (14.31 ) entails an implication 
An <2b ^ An+i < 2b and we can conclude that 

] 

b<6 \ ^An<26 WneN. (4.4) 
Ao<2b J 

Hence if we fix < e < ^ and assume that HmoIIw"! + I |^i-i/p^p^ + | I'W'ml |wpi(ri„) < e then 
starting the iteration from w'^) = (0, 0) we have 

+ \\w''\\w^ < 25 VnGN. □ (4.5) 

The next lemma almost completes the proof of the Cauchy condition in x L^(L2) for the 
sequence of iterated solutions. 

Lemma 12. Let the assumptions of Lemma [77] hold. Then we have 
where M is the constant from ( 14.71) . 



(r)' 
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Proof. Subtracting (ll.8l) rr. from (ll.Sh r, we arrive at 

- w"^^) - /iA(M"+^ - + ^)Vdiv - 

div - + 9xi(w"+i - t^^+i) + (m" + Mo) ■ V(i(;"+^ - = 

= It;") - 10™) + (m" - M™) ■ Vw;™, 

n ■ 2/iD(u"+i - M^+i) ■ Ti + / - u"'+^) ■Ti\r = 0, 
n- - w'^+^)|r = 0, 

The estimate (|2.4I) applied to this system yields 

||F(m",w'^) -F(u™,w'^)||y* + ||G(m",u;") -G(m'",u)'")||l2 + - u™) ■ Vio^Hl^- 

In order to derive (14.61) from the above inequality we have to examine the l.h.s. The part with G 
is the most straighforward and we have 

||G(«^^l;") - G{u"',wn\\L, < E{M) {\\u^ - u"^\\hi + |k" - ^"^IU^(l.)). (4.7) 

The function F is more complicated and we have to look at the difference more carefully. A 
direct calculation yields F(u", w"") - Fiu"^, w"^) = F"'"" + F^'"", where 

\\F^n\v* < E{M) (||«" - + ||«;" - ^/;"|U^(L.)) (4.8) 

and 

^2"'™ = -[(57r'(«;") - 57r'(«;'")]Vu;" + 57r'(ty'")V(u;" - w;™) =: F^f + Fa^f , (4.9) 

where 57r'(-) is defined in (|2.2I) . Since we are interested in the y*-norm of F^''", we have to 
multiply F21'' and F22'' hy v eV and integrate. With F2'^ we get 

(57r'(u)'") V(«;" - w™) ■ rfx = 
- 57r'(w'")(«;" - w™) div fix - /^(w" - w"") V [(57r'(w'")] ■ 

and thus we have to estimate 57r'(ty'") in terms of lo™. Using (|2.3I) we can write 

Sn'{w"') = w"' [ 7i" [tw"^ + l]dt, (4.10) 

what yields ||^7r'(i(;'^)||i^ < C(7r) || if"- 1| 2,^. Now we have to estimate ||V57r'(i(;'^)||ip. Since 
TT is a - function (and this is the only point where - regularity is needed) we can take the 
gradient of (14.101) and verify that \\^S'n:'{w"^)\\Lp < C{'n:)\\Vw"'\\Lp- Thus we have 

\J^57r'{w"'){w''-w"')divvdx\ < \\67i'{w"')\\l^ \\w'' - w'^Wl^ lldivi;!!^^ < 

<C\\w^\\w^ \\w^ - w^Wl^^l,) \\v\\v. ^ ^ 
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Next, since p > 3, by the Sobolev imbedding theorem we have 

\ J^iw"" -w'^)W[6n'{w'^)]-vdx\ < 
< \\w^ - w"'\\L,\\Vd7r'iw"')\\Lj\v\\L, <C\\w''^\\w^\\w^-w^\\l^(^l,)\\v\\v ^ ^ 

Combining (|4.11l) and (14.121) we get 

m^rWv* < E{M)\\w- -w^lL^^L^y (4.13) 
In order to estimate F2'^ we will use again fl2.3l) to write 

Sn'iw'') -dn^w"") = {w'' -w"^) [ p" [t w'' + {1 - t) w"" + 1] dt , (4.14) 

Jo 

what yields \\S7r'(w") — 67r'(w'^)\\L2 < C \ \w'^ — w'^\\l2- With this observation we can estimate 

\ Q6n'{w"') - S7i'{w"')]Vw'' ■ V dx\ < ||57r'(w") - 57r'(w™)|UJ|Vw;"||Lj|^^||L6 < 

what yields 

\\F^;r\\v* < EiM)\\w'' - w'^Wl^^l,)- (4.15) 
Combining the estimates on f]^''^,F^^^ and F^;^ we get 

IIFK,^/;-) - FK, w'")||v* < F(M) [ll^z" - w^Wh^ + \\w^- w'^Wl^^l,)]- (4.16) 

The part that remains to estimate is (u" — n"*) ■ Vto'". We shall notice here that this is the term 
which makes it impossible to show the convergence in x directly. Namely, if we would 
like to apply the estimate (12.351) to the system for the difference then we would have to estimate 
1 1 (m" — li™) • Viy™! Iwpi what can not be done as we do not have any knowledge about | |vi/2. 

Fortunately we only need to estimate the L2-norm of this awkward term, what is straightfor- 
ward. Namely, we have 

- u"^) ■ Ww'^Wl, < Ik" - u^\K\\Vw^\\l, < C\\w^\\w^ \\u^ - u^Whi, (4.17) 

since q = < 6 for p < 3. We have thus completed the proof of (14.61) . □ 

Now, lemma [TTI implies that the constant E(M) < 1 provided that the data is small enough 
and the starting point (u^,w^) = (0,0). It completes the proof of the Cauchy condition in 

X Loo(-^>2) for the sequence {(it", to"")}. 

Remark. Lemmas [TT] and [T2l hold for any starting point {u^,w^) small enough in x W^, 
not necessarily (0, 0), but we can start the iteration from (0, 0) without loss of generality. 
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5 Proof of Theorem [T] 

In this section we prove our main result, Theorem \T\ First we show existence of the solution 
passing to the limit with the sequence (m", ly") and next we show that this solution is unique in 



the class of solutions satisfying (11.31) . 

Existence of the solution. Since we have the Cauchy condition on the sequence (m", w") 
only in the space H^(Q)xLooiL2), first we have to show the convergence in the weak formulation 
of the problem (11.61) . transfering the derivatives of the density on the test function. The sequence 
(m", w") satisfies in particular the following weak formulation of (11.81) 



Jq{v ■ m""*"^ + 2/iD(M""'"-^) -.W V + z/ div w""*"-^ divv — 7 w"+^div v} dx ^ ^ 
+ Ir fi^"^^^ ■ '^i) ■ ^i) = Fiu"", ■vdx + /p Bi{v ■ n) da ^ ^ 

and 

- ^ [m" ■ V0 + div w" 0] dx + J^^^^ (j)da = 

= ^(^(m", u;") - div dx + Win (p da *^ 

V(^;, (p) eV X C°°{n), where = [1 + + uoY^\ + UoY^\ («" + UqY^^]. 

Now using the convergence in x Loo{L2) combined with the bound (14.11 ) in W.^ x 
we can pass to the limit in (15.11 ) - (15.21) . The convergence in all the terms on the r.h.s. of (15.11 ) is 
obvious and the only nontrivial step to show the convergence of F(m", w"^) is to show that 



/ 

Jn 



5'k'{w'^)Vw'^ ■ vdx / 5n'{w)Vw ■ V dx. 

Jq 



In 

To show the above convergence it is enough to verify that 



K(57r'(t/;") - wd7r\w)] div v dx ^ (5.3) 



n 



and 



[ KV57r'(ti;") - wV67r'{w)] -vdx^O. (5.4) 
Jq 

Applying again (12.31) we have 6tt'{w^) — 6tt'{w) = {w^ — w) 7r"{l+tw^ + (1 ~ i)^) dt, hence 

\\5'k\w'') - 5-n\w)\\L^ <C\\w''-w\\l^, (5.5) 
what implies directly (|5.3I) . To show (|5.4I) we integrate by parts arriving at 



«;"V57r'(«;") - wV5Tr'{w)] ■ v dx = - / {dn'{w'') - 6ti'{w)) divvdx 
n Jq 

- I {5ti'{w'') - 5n'{w)) V ■ Vw" dx+ f (u)" - w) V5tt'{w) ■ v dx. 
Jq Jq 

Now all the terms on the l.h.s. converge by (15.51) . hence we have shown (|5.4I) . 
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We conclude that (m, w) satisfies 

Jq{v ■ dx^u + 2fi'D{u) : V V + z/ div u divv — •yw div v} dx ,^ 
+ jp f{u ■ Ti) {v ■ Ti) da = F{u, w) ■vdx + /p Bi{v ■ n) da 

V f G y . In (15.2p we have to check the convergence in the boundary term. We can use the same 
argument as in the proof of the existence of solution to the linear system when we have passed 
to the limit with finite dimensional approximations. Namely, in fact ly" satisfies the Cauchy 
condition not only in Loo{L2). A stronger fact holds that w"^ is a Cauchy sequence in L2{Vlx^) 
for every xi G [0, L], where Vt^^ denotes the Xi-cut of Vt. In particular ui" — > C in ^2 (rout) for 
some C G L2{Tout) and since sup^^g[o,i,]ll'"^lli2(n,,) < 00 we conclude that C = w|ro„f This 
result combined with the obvious convergence of other terms in (15.21) implies 



— / w[u ■ V (j) + div ucj)] dx + I w(f)da= I (j){G{u, w) — div u) dx + / Win(f)da (5.7) 

V0 G C^in), where u=[l + {u + uq)^^\ {u + uo)^'^\ {u + uo)^% 

Hence we have shown that (tt, w) satisfies (15.61) - (15.71) . the weak formulation of (11.61) . Now 
we want to show that the strong formulation also holds. 

The bound in x implies (-u"'=,w"'=) ^ {u,w) in x for some {u,w) G 

X Wp. On the other hand, we have {u"''' ,w^'=) — > (u, w) in x Loo{L2), thus we conclude 
that (u, ti?) = (tt, If). 

Hence we can integrate by parts in (15.61 ) - (15.71 ) to obtain 

/f7 [^(^^ ~ l^^u — + z/) Vdiv u + 'jVw] ■ V dx ^ ^ 

= Jp [Bi{v ■ Ti) — n ■ [2/iD(u) + z/divitld] ■ v — f{u ■ Ti){v ■ Tj)] da 



and 

/ [wx-i^ + {u + uq) ■ Vw](p dx = / [^(Ujiy) — divM]0(iz. (5.9) 

From these equations we conclude that (ll.6l) i 9 are satisfied a.e. in f2 and i\l.6\i -^ is satisfied a.e. 
on r. It remains to verify that (II. 6L is satisfied a.e. on T and (11.61) .^ holds a.e. on Fj^. The 
condition (|1.6I )zl results from the convergence — > m in if^. 

Finally, ^ w in implies that w^"|r,„ ^ ^'^'"^|ri„ in ^p(rjn)- On the other hand 
^"^bin '<J^in in (Fjn) since it is a constant sequence. We conclude that w\Yi^ = Win- 

Uniqueness. In order to prove the uniqueness of the solution consider (vi, pi) and (t'2, P2) 
being two solutions to (11.11) satisfying (11.31 ). We will prove that 

11^1-^^211^1 + IIpi -P2IIL =0- (5-10) 
For simplicity let us denote u := vi — V2 and w := pi — p2- We will show that 

< -ElltyllL, (5.11) 
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and 

IkllL. <C||^/||^,i, (5.12) 
what obviously implies (|5.10l) . Subtracting the equations (11.11) for (fi, pi) and {v2, P2) we get 

w V2 ■ Vv2 + piu ■ Vv2 + PiVi ■ Vu — fiAu — (/i + z/) Vdiv u + It,Vw + wVI-,, = 0, 
Pi div u + w div V2 + u ■ Vp2 + Vi ■ Vw = 0, 
n ■ 2p'D{u) ■ r|r = 0, 
n ■ u\r = 0, 

Mr.n = 0, 

(5.13) 

where ^ 

[ 7v'{{tpi) + il-t)p2)dt. (5.14) 

Notice that I,, E Wp since pi E Wp and n E C^. In order to show (15.1 II) we follow the proof of 
(12.41) multiplying (I5.13l) i by pi u (it will be clarified soon why take the test function pi u instead 
of u). Using (12.61) we get 

j {2p'D^{u)+upidiv'^u)dx+ J ^2p[{pi - 1)D(m) : Vw + D(w) : (m (g) Vpi)] + u {div u)u ■ \/ pi^ dx 

h 

— J I M Vpi + p^M^ ■ Vv2 + uw P1V2 ■ Vf 2 1 (ix + J p\ {vi ■ Vu) ■ u dx 



+ / piw u -V In dx — I wu ■ \/{I.n-pi) dx — / It,w Pi div u dx + / pifu'^da = 0. 
in Jq Jn Jr 



/4 /5 



We have + I/2I < £^(|| w| 1^1 + 1 1 1 1^2) ^^'^ order to deal with let us split it into two 
parts: 

2/3= / {{plv'i'' -l)d^^\u\'^ + pIv'^^ d^M^ + pIv^{^^ d^,\u\^}dx+ [ d^,\u\'^ dx . 
Jq Jn 



We have I/3I < and /| = J^\u\'^n^^^ da = — /p {ul"^ da + J^^^Jul"^ da. In order to 

examine I4 and we have to differentiate (15.141) what yields 

V4 = J^Vpi + />P2, (5.15) 

where 



Il= [ Tx"{tpi + {l-t)p2)tdt and Il= f Tx"{tpi 
Jo Jo 



+ {l-t)p2)il-t)dt. 
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We have 



/ pilluw\/pidx\ < WpiIIWl^ llVpilli^ \\u\\Le \ \w\\l2 < EiWuWjj, + 
Jn 



\w\ 



and the same for pi I^uwV p2 dx. Thus the application of (15.151) to L yields \Ia \ < E {\\u\\'jji + 
1 1 1^1 11^). To estimate I/5I it is enough to use (15.151) to compute V^Ij^pi) and then with the same 
arguments as in case of I4 we get I/5I < E + HwHl^)- Summarizing our estimates we 

can write 



\u\ 



< It, w pidiv u dx + E WwWj^^. 



The boundary integrals over Tq and Tout will be nonnegative for any / > and the integral over 
Tin will be nonnegative for / large enough on Fj^. Now in order to obtain (|5.1 II) from (15.161) we 
can express pi div u in terms of w using the equation (15.131) 9 (this is why we have tested (|5.13|) i 
with pi u instead of u) and rewrite (15.161) as 

||w||^i < — / InW^divv2dx— / lT,wu-Vp2dx— / It,viw ■ Vw dx +E\\w\\^i^^. (5.17) 

Jn Jn 



h I7 h 



We verify easily that l/gl + I/7I < -E" (| |m| 1//^ + | l^^l lij)- We have to put a little more effort to find 
a bound on Jg. Let us integrate by parts: 

2/3= / Ij.Vi'Vw'^ dx = — / w'^div{lT,vi) dx + / w'^It.Vi ■ nda. 
Jn Jn Jr 

The boundary term reduces to Jp ^ I-^w'^v^^ da > and in order to deal with the first term on 
the l.h.s. notice that 

div(/^ vi) = divvi h + ll^i- Vpi + vi ■ Vp2, 

hence 

2/8 < — / nP'divvi I.„dx— / w"^ Vi ■ Vpi I^dx — / w"^ Vi ■ Vp2 dx . 
Jn Jn Jn 



Tl t2 j3 

-"8 -"8 



Obviously we have < i? | |w| H^. In order to bound /| we can apply the continuity equation 
that yields i/j ■ Vpi = —p^ div fj, what implies |/|| = \ J^w'^ pidivvi I^dx] < E \\w\\1^.lnt\ie 
term /| we can rewrite the mixed component as vi ■ Vp2 = m ■ Vp2 + '^^2 ■ Vp2 and conclude that 
\Ii\ < -E'dkll^i + ll'U'lliJ- Combining the above results with (15.161) we get (15.111) . 
In order to show (15.1 II) we express the pointwise value of w using (15.131) 9: 



w'^(Xi,X2) 



/ wWs{s,X2)ds = — / —tytW div u{s, X2)ds 
Jo Jo f I 
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Jo v\ ^ ' ^ Jo v\' 



2 I 2 I 2 



We estimate directly the first two components of the l.h.s. obtaining 

\ w\dx <t\\w\f^^^C{e)\\u\^Hr Ve>0 

and ^^w\dx < E (WwWj^^ + | |m| |^i)- To complete the proof we have to find a bound on w^. To 
this end notice that 



w- 



dx = - f I — rp- \vf^dx^w'^ + vf''dx^w'^~\ dxdxi, 
2 Jo Jp^^ v[ 



where P^^ = VIq x (0, xi). Integrating by parts in the inner integral we get 

(2) ,,(3) 



n 



dx=U {-/ w'[d^,\r + d^,\]dx+ ^[v^M'^+v['^n^'^]da}dx^. 
2 Jo Jp,,-^ vl' V{' JdP^-^ vl ' ^ 



\w\ 



The boundary integral reduces to Jy^^^qp w'^ v ■ nda = 0, what implies J^w^dx < E\ 
and (15.121) easily follows completing the proof of the uniqueness, and hence the proof of the 
Theorem. □ 



6 Appendix 

Vorticity on the boundary. In order to show the boundary relation (|2.1 Ih /i we have to dif- 
ferentiate (11.91) 4 with respect to tangential directions at a given point xq E T. Without loss of 
generality we can assume that n{xo) = (1, 0, 0), ti{xo) = (0, 1, 0) and T2{xo) = (0, 0, 1). Then 
we can rewrite i\l.9h as (all the quantities are taken at xq): 

J /i(M^2+^^^l) + = 5i, .... 

\ fi{u\,+u^,) + fu' = B2. ^""-'^ 
Differantiating (1 1.91 ) a with respect to the tangential direction ti we get 



(—n)-u + u\2=0. (6.2) 



If we denote by xi the curvature of the curve generated by ri then we have = xi^i ^^'^ 
(16.21) can be rewritten as Xi(ri ■ u) + u^,2 = 0. Combining this equation with (I6.1l) i we get 



u^a -u\2 = (2xi ){u- n) + — , 

/i /i 

what is exactly (12.1 Ih . (12.1 ll) ^ can be shown in the same way differentiating (11.91) 4 with respect 
to the tangential direction T2. 
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Lemma 13. (interpolation inequality): 

Ve > 3C(e, fi) such that V/ G W^p^Q); 

||/||L,<e||V/||i, + C||/||z.,. (6.3) 

Proof. Inequality Q results from the inequality 1 1/| < C{p, ^]) 1 1/| 1 1/| for 2 < 

j9 < oo, where = (see fT], Theorem 5.8). Using Cauchy inequality with e we get 16.31 □ 

The last auxiliary result we use is a following fact on finitely dimensional Hilbert spaces (the 
proof can be found in [|24l ): 

Lemma 14. Let X be a finite dimensional Hilbert space and let P : X X be a continuous 
operator satisfying 

3M>0: (P(O,O>0 for Ml = M. (6.4) 

Then^C ■■ 11^11 < ^ and P(^*) = 0. 
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